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Helical and helicoidal precipitation patterns emerging in the wake of reaction-diffusion fronts 
are studied. In our experiments, these chiral structures arise with well defined probabilities Ph 
controlled by conditions such as e.g., the initial concentration of the reagents. We develop a model 
which describes the observed experimental trends. The results suggest that Ph is determined by 
a delicate interplay among the time- and length-scales related to the front and to the unstable 
precipitation modes and, furthermore, the noise amplitude also plays a quantifiable role. 
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Helices and hclicoids are present from nano- to macro- 
scale (ZnO nanohelices [l|, macromolecules and inorganic 
crystals with a helical structure 0, Q , precipitation he- 
lices fiber geometry of heart walls [7j). Forma- 
tion of these fascinating structures generally follows two 
routes. First, templates with chiral symmetry (e.g., ora- 
gogel fibers) may exist in the system, and the symmetry 
is just transcribed to a structure (e.g., inorganic materi- 
als Q) at a larger scale. Second, spontaneous symmetry 
breaking may occur through the self-assembly of achiral 
building blocks into a helical/helicoidal structure, as e.g. 
in case of crystals with chiral morphology [D, 0| • 

Theoretically, the symmetry-breaking route is more in- 
teresting. Universal aspects may emerge and the robust 
features of this self-organization process may be impor- 
tant for applications as well. Indeed, control over creat- 
ing helical structures would make engineering (in partic- 
ular, the bottom- up design of micro- patterns |l0() more 
flexible since chiral morphology of materials are known 
to affect their physical (electronic) properties [H [ll[ . 

In order to develop insight into the genesis of he- 
lices/helicoids, we investigate an emblematic example 
of pattern formation, namely the formation of precipi- 
tation patterns in the wake of reaction-diffusion fronts 
2 13 1- The motivation for this choice comes from the 
observation that helicoidal structures have an axis, and 
the correlations are simple in the plane perpendicular 
to the axis. This suggests that building the perpen- 
dicular correlations in the wake of an advancing pla- 
nar front may be a simple and natural mechanism of 
creating heliccs/helicoids. Additional motivation comes 
from the existence of a large body of knowledge in the 
related Liesegang phenomena [HI El- It allows the 
use of well-established experimental and theoretical ap- 
proaches, thus making it easier to develop a novel view 
on the formation of helical structures. 

Our main results concern the probabilistic aspects of 



the symmetry-breaking route. We determine the proba- 
bility Ph of the emergence of single helices/helicoids in 
Liesegang-typc experiments as the conditions such as the 
initial concentration of inner or outer electrolytes, or the 
temperature are changed. Ph is found to be well repro- 
ducible and large (Ph > 0.5 for some parameter range). 
The results are understood by expanding and simulat- 
ing a model of formation of precipitation patterns fl3 |. 
We explicitly observe that the origin of helices/helicoids 
is not to be found in the fluctuations and asymmetry of 
the initial- or boundary conditions [HI, [n| . Instead, the 
growth of unstable modes, the dynamics of the front, and 
the bulk fluctuations (noise) combine to yield the helices. 
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FIG. 1: Regular Liesegang- (leftmost tube) and helicoidal pat- 
terns (all other tubes) in agarose gel with the numbers corre- 
sponding to the tube radius R measured in mm. R is varied 
at fixed experimental conditions (T=22°C, [Cu 2+ ]o=ao=0.5 
M, and [CrO4~] =&0=0.01 M) and the probability of helicoid 
formation Ph is displayed (right panel). No helicoid appears 
for R< R c . 

In our experiments, we study the precipitation reaction 
Cu 2+ (aq) + Cr04 - (aq) -> CuCr0 4 (s) in 1% agarose gel. 
The gel soaked with K 2 Cr04 (inner electrolyte) is placed 
in a test tube and a solution of CuCl2 (outer electrolyte) 
is poured on top of the gel. Setting the concentration of 
the outer electrolyte an order of magnitude larger than 
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that of the inner electrolyte yields a reaction front diffus- 
ing into the gel, and a Liesegang pattern of precipitation 
bands forms behind the front (FigfT]). Frequently, how- 
ever, helicoids evolve from the same macroscopic exper- 
imental conditions (FigfT]). We quantified the stochastic 
nature of this intriguing phenomenon by varying the con- 
centration of the outer (ao) and inner (bo) electrolytes, 
the radius of the test tube (R), and the temperature (T), 
and measuring Ph using 10 independent experiments for 
each parameter set. 

Similar experiments were carried out in a quasi two- 
dimensional geometry as well. The gel (with the inner 
electrolyte B) was placed in the gap between two test 
tubes of slightly different radii (5R = 2 mm), thus effec- 
tively confining the pattern to the surface of a cylinder 
(FigfSJ). In this geometry, we observed the formation of 
regular Liesegang rings, single helices, double helices, and 
more complex patterns for large R. 
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FIG. 2: Transforming the thin layer of gel in the tube-in- 
tube experiment into a two-dimensional strip. The Liesegang 
bands and the helices in the experiments were obtained in 
agarose gel at T = 22 °C, a = 0.5 M, and b = 0.01 M, with 
the radii of the outer and inner tubes being 8 and 6 mm. 
The scaled parameters used for the simulations (columns on 
the right with the precipitate shown in white) were ao = 80, 



bo = 1, o = 0.8, A = 0.2, r? = 0.05, and L y 



64. 



FigfT] shows regular bands and helicoidal patterns in 
test tubes of various radii, together with the measured 
probability of helicoid formation (Ph)- We observe no 
helicoids below a critical radius (i? c =1.5 mm), in agree- 
ment with theoretical expectations based on a simplified 
model where the reaction front moves with fixed velocity 
151 Il6j . For R > R c , one finds that Ph increases with 
increasing R, and reaches rather large values (Ph ~ 0.7 
at R = 8 mm) before decreasing again. The decrease 
at large R is due to the noticeable proliferation of com- 
plex structures (double or triple helicoids, disordered pat- 
terns) which suppress the weight of single helicoids. 

Before describing the experiments further, we turn to 
the theory since it allows a more concise discussion of 
the results. Theories of Liesegang-type patterns combine 
the properties of a moving front (i.e. where and at what 
rate the reaction product, A + B — s- C, appears) with the 



details of the precipitation (i.e. how the reaction prod- 
uct, C, turns into precipitate). While the front proper- 
ties have been thoroughly studied and understood both 
theoretically an d experimentally 19 L 20 L the dy- 



namics of precipitation is more debated ImTiI- The 
competin g p re- and post-nucleation views can be com- 
bined [H , 21 1 , and we shall use a simple version [13 based 
on the Cahn-Hilliard equation with noise added |22l - E3 | . 
This equation features spinodal-decomposition-type fast 
dynamics, as well as slower, nucleation-and-growth pro- 
Driving it with a reaction zone gives us a flex- 
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cesses 

ible model with a variety of pattern-formation regimes. 

The reaction front appears due to a strongly inhomo- 
geneous initial distribution of the reagents A and B. 
The reaction takes place in a gel (occupying the half 
space x > 0) and, initially, the inner electrolyte B is 
homogeneously distributed [b(x > 0,y,t = 0) = bo]- 
The outer electrolyte A of much higher concentration 
[a(x < 0, y,t = 0) = ao with ao 3> bo] is brought into 
contact with the gel at t = 0. Assuming a second-order 
irreversible reaction A + B — > C, the front invading the 
gel can be described by the equations 



dta = P^Aa — kab 
dtb = D B Ab-kab 



(1) 
(2) 



where both the reaction rate k and the diffusion coeffi- 
cients, which are assumed to be equal (P^ = Dg = D), 
are set to 1 by an appropriate choice of the time- and 
length-scales [26[. The front, specified in terms of the 
rate of production of Cs (kab), is narrow and moves into 
the gel diffusively (the position is given by Xf — yj2D~fi 
where Df is a function of P and &o/ a o)- The front 
leaves behind a constant concentration of Cs (co), where 
Co depends on P and 6o/ao, and it is practically inde- 
pendent of k. Provided the system with Co is unstable 
or metastable, a phase separation of Cs into regions of 
high (ch) and low concentrations (q) takes place. This 
process is described by the Cahn-Hilliard equation with 
source (kab) and noise (rj c ) terms added 



dtm = — AA(m — m 3 + a Am) + kab + r\ c 



(3) 



Here m is the concentration of Cs shifted by c — (c^ + 
ci)/2 and scaled by c = (c.% — Ct)/2, so that m = (c — c)/c 
is 1 for c = Ch and m = — 1 for c = eg. The parameters 
A and a are the rescaled kinetic coefficient and surface 
tension, respectively [3, 23, 24 1. Their ratio t u w ct/A 
defines a characteristic timcscalc of the growth of unsta- 
ble modes in precipitation. Comparing t u with the time 
the front passes through a region determines whether 
slow nucleation-and-growth or fast spinodal decomposi- 
tion dominates the pattern formation. 

Adding noise (r] c ) is essential since the formation of he- 
lices is a symmetry-breaking process. Furthermore, the 
noise widens the available regions of the meta- and un- 
stable states (see earlier morphological phase diagrams 
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FIG. 3: The probability of helicoid /helix formation Ph in experiments (large red dots) and in simulations (small symbols). 
Displayed are the dependence on the outer- and inner electrolyte concentrations (Panels a and c), on the noise amplitude (b), 
and on the radius (width) of the system (d). The values of the parameters kept fixed in a panel, and the experimental estimate 
of the scale factor c are discussed in the text. Statistical errors are shown for the experiments and for a single set of simulations. 



of Liesegang patterns [I?}). Noiseless Cahn-Hilliard type 
dynamics where the front moves with fixed velocity have 
been much studied [H,^. In these cases, however, noise 
was present in the initial state, and complex morpholo- 
gies resulted from complex initial conditions or from com- 
plex motion of the reaction front. Our model without the 
noise reproduces the properties of the regular Liesegang 
patterns [3, . Inclusion of the bulk noise allows us to 
demonstrate the existence of helices, and understand the 
experimental trends in their emergence. 

From a theoretical point of view, the tube-in-tube ex- 
periments are the easiest to describe. We can cut and 
open the cylinder as shown in FigJ5] and treat the thin 
layer as a two-dimensional strip of width L y = 2-7rP, and 
length equal to the tube length L x [3lj |. Accordingly, 
eqs.([l][3|) are solved in a rectangle of size L x x L y with pe- 
riodic boundary conditions in the y direction and no-flux 
boundary conditions at the lower edge (x = L x ,y). At 
the upper edge (the initial location of the front (x = 0, y), 
we use a slightly idealized boundary: the concentration 
of the outer electrolyte is kept at a constant value a(x = 
0,y,t) = ao/c while no- flux condition is adopted for B 
and C. The initial conditions reflect the experimental 
setup: b(x > 0, y, t = 0) = bo/c, a(x > 0,y,t = 0) = 0, 
and c(x,y,t = 0) = 0. The discretized noise term r\ c is 
implemented by moving Cs between neighboring sites at 
a rate rj c = y/cr where r is uniformly distributed in the 
interval [—77,77]. In the following, 77 is called the ampli- 
tude of the noise. 

Our simulations indicate that both the Liesegang 
bands and the helices emerge in a wide range of the pa- 
rameters. There are, of course, some constraints, e.g., 
77 must be sufficiently small for the phase separation to 
take place. Examples of simulations are shown in Figf5] 
(rightmost two columns), where a Liesegang pattern and 
a helix are displayed [26| . A general feature of the simula- 
tions is that the chirality of the helices is random within 
the statistical error of 100 independent simulations. This 
is in agreement with the experiments where, out of 96 
helicoids/helices, the ratio of left- and right-handed pat- 



terns is 50/46. We consider this as experimental evidence 
that the macroscopic symmetry breaking is not driven by 
microscopic objects of given chirality. 

To characterize the emergence of the helices quanti- 
tatively, we collected data by varying ao, bo, f] and L y , 
and determined Ph from the outcome of 100 simulations 
with distinct random number sequences for r\ c . Since the 
kinetic coefficients A and a cannot be controlled exter- 
nally, we kept them fixed (A = 0.2, a — 0.8) throughout 
the simulations. 

First, we varied ao/c and 77 while keeping bo/c = 1 and 
L y = 64 fixed. FigsOk,b shows that Ph is remarkably 
large, it increases with ao and reaches Ph ~ 0.4 — 0.6 
for large ao/c. Similar trend is also seen in the exper- 
iments. Since ao/bo determines the front motion, with 
larger ao corresponding to faster diffusion, an important 
conclusion from Figj3K is that fast motion of the front 
facilitates the emergence of helices. 

Figj3^.,b also show that no helices form even for larger 
ao /c if the noise is too small. Increasing the noise first 
increases Ph, then Ph saturates in the region 0.05 < 
77 < 0.09 and, finally, ft -> due to the absence of 
phase separation above 77 rj 0.09. Comparing FigJSjD with 
experiments is difficult since the link between 77 and T 
is through complex changes in diffusion, reaction rates 
etc. Our experiments indicate that Ph increases with T. 
This is in agreement with Figj3b provided 77 ~ T and the 
experimental T corresponds to small values of 77. 

We also varied bo/c and 77 while fixing ao/c — 80 
and L y = 64 (Fig[3j:). The probability Ph was found 
to be maximal in the middle of the spinodal region 
(bo/c w 0.9—1.1) where isotropic precipitation structures 
develop through fast-growing, linearly unstable modes. 
Comparing the simulations (Fig(3J;) with experiments is 
difficult since neither 77 nor the experimental concentra- 
tion scale, c = (ch + q)/2 rj c/j/2, are known. We esti- 
mated Ch by assuming that all the precipitate was in the 
helices and all the Bs reacted and turned into Cs. This 
estimate left an apparent shift between the experimental 
and simulation points (FigJSJ;). The shift may well be 



the consequence of overestimating ct (e.g., not all the Bs 
reacted, or the bands arc wider than their optical width). 

The effect of increasing width (L. y ) is displayed in 
FigOJi. The experimental parameters are described in 
Fig|TJ while in simulations, we used ao/c = 15, ba/c = 1, 
and rj = 0.04. The experimentally observed lower thresh- 
old for the emergence of helices is clearly present (L yc w 
32), and one can also recognize the trend that Pa first 
increases with L y and then decreases for large L y . As 
in experiments, Ph decreases due to the proliferation of 
more complex structures. Complexity builds up for large 
L y since more long-wavelength transverse modes (modes 
in the y-direction) can fit into the system. They are un- 
stable modes of the Cahn-Hilliard dynamics facilitating 
the formation of more intricate patterns. 

The common trends found in experiments and simula- 
tions suggest that our model contains the right ingredi- 
ents, and we can develop a physical picture of helix for- 
mation by observing the simulations. Fig |4] displays two 
examples of time evolutions with parameters set to have 
roughly equal probabilities for bands and helices. There 
are many ways of choosing such parameters but the char- 
acteristic features of the dynamics are always the same. 
Essential among them is that, initially, the reaction front 
moves fast enough to produce a domain where the sys- 
tem is unstable and roughly homogeneous (fuzzy regions 
in the t = 720 — 960 plates in FigfJ]). The homogeneity 
makes possible the generation of isotropic patterns which 
compete with the anisotropic influence of the front fa- 
voring band-formation perpendicular to the front motion 
(t = f440 - f 920 plates in Figg]). The outcome of this 
competition determines whether Licsegang bands, single 
helix, double helix or more complicated patterns emerge. 

One can quantify the above picture by noting that ho- 
mogeneous patterns can form only if the front moves a 
distance of the order L y in a time, 77 = L y /2Df that 
is smaller than the time, r„, required for the precip- 
itation to develop. To estimate t u , we calculate the 
growth rate, cj&» = A/4cr rj 1/t u , of the fastest grow- 
ing mode of wave- number k* = Y/\/2a using linearized 
Cahn-Hilliard dynamics for a quench to the middle of 
the miscibility gap [ro(0) ~ 0]. Then, assuming that 
the homogeneous structure emerges from the noise, we 
have y/rjexp (u>k*T u ) m(r„) sa 1, and the inequality 
Tf < t u yields an upper limit for the width of a sys- 
tem L y < 4Dftj\ lnryl/A where helix can form. A lower 
limit can also be found since the characteristic size of the 
domains (L* ps 2-rt/k* = 27tv / 2ct) formed by the fastest 
growing modes must be smaller that the width of the 
system {L* < L y ), otherwise no structure forms in the y 
direction. The combination of the two inequalities 

8ir 2 a < L\ < 4D f a\\nr)\/A (4) 

reflects some of the trends observed in the experiments 
and simulations. Namely, the formation of helices are 
facilitated by a fast moving front, i.e. by Df being large 
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FIG. 4: Time evolution of the precipitate (white regions) for 
parameters ao/c = 60, bo/c — 1, 7/ = 0.02, and L y — 64. 

which, in turn, requires ao to be large and, furthermore, 
there is a minimal width below which no helices form. 

Since the width {L y ) is bounded from both sides, it 
may happen that no helices can form. When searching 
for helices one should, in general, use a fast front (e.g., by 
selecting large ao) and create an unstable state behind 
the front by placing the system deep in the miscibility 
gap (e.g., by experimenting with bo). Finding the right 
temperatures is also important but it is a rather complex 
problem left for future studies. 
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SUPPLEMENTARY INFORMATION 



A. Detailed model description 

Liesegang patterns are formed in the wake of moving 
reaction-diffusion fronts. The reaction front emerges due 
to the inhomogcneous initial distribution of the reagents 
A and B. Namely, the reaction takes place in a gel occu- 
pying the half space x > where, initially, the inner elec- 
trolyte B is distributed homogeneously [b(x,y,t = 0) = 
bo6(x)]. The outer electrolyte A of much higher initial 
concentration [a(x, y,t = 0) = ciq9(— x) with clq ^> bo] is 
brought into contact with the gel at t = 0. Assuming a 
second-order, irreversible reaction A + B — > C, the front 
invading the gel is described by the equations 



dta = DyiAa — kab 
d t b = D B Ab-kab 



(5) 
(6) 



where k is the reaction rate and, for simplicity, the dif- 
fusion constants of the reagents are assumed to be equal 
(Da = Db = D). This assumption will be used through- 
out the two-dimensional simulations. 

The front, characterized by the spatio-temporal prop- 
erties of kab (rate of production of Cs), has been stud- 
ied in detail [It], EH ■ It is narrow and moves diffusively 
(its position is given by Xf = y/2Dft where Df can be 
expressed through D and bo/ao). It leaves behind a con- 
stant concentration, c = Co, of the reaction product C, 
where the parameter cq is determined by D and bo/ao, 
and is practically independent of k. 

Assuming that no intermediate complexes are formed, 
the next stage of the formation of the precipitation pat- 
tern is the separation of the reaction product, C, into 
high- and low-concentration phases. At a coarse grained 
level, the phase separation can be d escribe d by the gen- 
eralized Cahn-Hilliard equation 



14, 22-2 



d t c = X A(Sf/Sc) + kab + rj cQ . 



(7) 



Here Ao is a kinetic coefficient, / is the free energy driv- 
ing the phase separation, kab describes the creation of 
C particles by the front, and rj C Q represents noise effects 
(thermal fluctuations, inhomogcncities in the gel, etc.) 
which conserve the total number of C particles. 

In order to describe the phase separation, the free en- 
ergy /(c) should have two minima corresponding to the 
low- (cf) and high (ch) concentrations of Cs in homoge- 
neous equilibrium states. It should also have a surface 
tension term preventing the formation of singularities at 
interfaces. As a convenient form with minimal number 
of parameters, one can take / as a Landau-Ginzburg free 
energy which is symmetric about 5= [ch + q)/2 



f(c) = --(c- 



c) 2 + l(c-c)* + ^(Vc) 



(TO 



where e, 7, and ao are phenomcnological parameters, 
and the minima of /(c) are fixed at Ch and ct by setting 
y/s/'y = (ch ~ Cf)j1 w where we use the fact that 
Ch ^> Q i.e. the gaps between the bands have very low 
steady-state concentration of Cs in the usual Liesegand 
experiments. 

Measuring concentration, time, and length in units of 



Ch - C£ 



Ch 

2 



(9) 



and, furthermore, making a shift in the concentration of 
Cs 



c- (c h + c e )/2 ^ c 
c 



(c h - ct)/2 
one obtains a simple set of equations 



(10) 



d t a = Aa — ab 
dtb = Ab-ab 
dtm = — XA(m - 



(11) 
(12) 

a Am) + ab + r) c , (13) 



(8) 



where A = Xos/D, a = aokc/De, r\ c = rj c o/kc 2 are 
the rescaled kinetic coefficient, surface tension, and con- 
served noise, respectively. 

A few comments are in order about the random as- 
pects of the dynamics. First, we note that randomness 
is not added to the reaction-diffusion equations (|11I12|) 
since the noise has been shown to be irrelevant in de- 
termining the properties of the A + B — > C type re- 
action fronts in the physically relevant dimensions [32| . 
Second, we recall that the noise term r] c in the gener- 
alized Cahn-Hilliard equation (Model B 23() describes 
the local concentration-fluctuations resulting from diffu- 
sive random motion of Cs. This noise conserves the total 
number of particles and is expected to be proportional 
to •Jc. Third, we point out that the amplitude of r/ c in 
near-equilibrium dynamics is proportional to the temper- 
ature and it is related to the kinetic coefficient A. How- 
ever, we have here a far from equilibrium situation, and 
no fluctuation-dissipation relation connects these param- 
eters. Accordingly, we shall treat A and the amplitude 77 
of r\ c as independent parameters. Of course, one expects 
that temperature is related to the noise and, in general, 
77 increases with temperature. In the simulations, the 
discretized noise term \r\ c in (Tl3|) ] was implemented by 
moving Cs between neighboring sites at a rate propor- 
tional to rjc = y/cr where r is a uniformly distributed 
random number from the interval [—77,77], and rj is the 
parameter which is called the amplitude of the noise. 

The above formulation is remarkable in that the num- 
ber of parameters (A , a , 77) is small compared to that 
found in the usual models of Liesegang phenomena. 
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Among the parameters, a does not appear to be im- 
portant since it just determines the width of the tran- 
sition between the high- and low-concentration regions. 
The parameter A, on the other hand, does play an es- 
sential role since it sets the timescale of the precipita- 
tion processes. Comparing this timescale with the time 
the front passes through a region determines whether the 
nucleation-and-growth or the unstable growth (spinodal 
decomposition) modes dominate the phase separation dy- 
namics. Finally, the noise r\ c is also important. First, 
because the formation of helices is a symmetry breaking 
process which does not happen without the presence of 
noise. Second, because it determines the phase diagram 
(the meta- and unstable regions) for the given parameter 
values of the system. 

Unfortunately, the parameters A ,a ,rj are not easily 
controlled in experiments. For example, when trying to 
amplify rj by increasing the temperature, one immedi- 
ately realizes that there are a number of parameters (dif- 
fusion coefficients, reaction rate, etc.) which have strong 
but largely unknown temperature-dependence with un- 
predictable combined effects. 

The parameters which are contollable in the experi- 
ments come from the initial preparation of the systems. 
They are the initial concentrations of the electrolytes 
(do , bo) and, furthermore, it turns out that the radius R 
of the test tube also sets some constraints on the emer- 
gence of helices. 



B. Additional information about the simulations 

The solution of the discretized equations (|1HI13P with 
the above boundary and initial conditions were carried 
out using a uniform grid with various combination of 
scaled L x and L y from the ranges of 32 < L y /l < 512 
and 512 < L x /l < 2048. The equations were integrated 
in time by the simple Euler method (fast and extensive 
search in the parameter space was feasible by using the 
parallel programming possibilities of video cards). For 
the results quoted and displayed in the paper, the grid 



spacing and the time step were 1.0 and 0.02, respectively. 



C. Three-dimensional simulations 

Three-dimensional systems can also be studied by sim- 
ulating the d — 3 versions of equations ([5][7J) . There are 
two changes compared to the d = 2 case. First, in order 
to be in agreement with the experimental setup, the peri- 
odic boundary conditions are replaced by no-flux bound- 
ary conditions in the direction perpendicular to the direc- 
tion of motion of the front. Second, we lift the restriction 
of Da = Db in our search for helicoids, and this means 
that an additional parameter 8 — Da/Db appears in the 
scaled equations. 

In the simulations, we find that both the Liesegang and 
the helicoidal patterns observed in the experiments can 
be qualitatively reproduced (see Fig{5]) . Unfortunately, 
the time-scale of the simulations compared to the d = 2 
ones is multiplied roughly by 2R where R is the scaled 
radius of the tube (the grid spacing used was 1.0). As a 
consequence, the computer power presently available to 
us is insufficient for obtaining good quality statistics for 
the probability of hclicoid formation in d = 3. 



i 


i 











FIG. 5: Precipitation patterns (Liesegang on the left, he- 
licoidal on the right) obtained in d = 3 simulations using 
the following scaled parameters ao = 100, bo = 1, a = 0.8, 
A = 0.2, t) = 0.02, 9 = 1.2, and R = 24. The only difference 
in the simulations is the initial seed for the random number 
generator. 



